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Abstract 

Channeling of recoil ions in various crystals strongly depends on their direction. This, along 
with the WIMP wind due to the rotation of the galactic disc generates daily modulations of the 
channeling. Since channeling affect the characteristics of the WIMP induced recoil events these 
modulations may be observed in single crystal detectors and serve as a WIMP signature. We 
suggest in particular searching for such modulations in the DAMA scintillation data. 

PACS numbers: 
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I. Introduction 



Astrophysical evidence strongly suggests non-baryonic matter constituting most of the 
matter's total cosmological mass density. Direct searches for such dark matter particles 
in low background, large and very sensitive underground detectors has been ongoing for 
several decades. Signatures separating signals from background can help this search via rare 
nuclear collisions of the Weakly Interacting Massive dark matter Particles (WIMP) in these 
detectors. 

One such signature, a ~ 10% annual modulations was suggested some time ago in a 
beautiful paper (DFS).jlJ] DFS assumed isotropic WIMP velocity distribution in the halo 
frame. Since the solar system is rotating with the galactic disc at a velocity W of 230 Km/sec, 
WIMPs impinging on earth have this net wind of velocity W. On June 2nd a component 
comprising ~ 1/2 of earth's rotation velocity of 30 Km/sec adds up to W increasing the 
WIMP flux, energy and counting rate and the opposite happens around December 2nd. 
Such modulations have been seen in DAMA's Na-I scintillation experiment^ but not in 
other WIMP searches 

"Channeling" , namely, the tendency of atoms, ions or electrons entering a crystal to move 
in channels along crystal axes or Bragg planes is well known and studied. See, eg., Refs. 



and 



io|. 



Channeled recoils collide less frequently and deposit more energy in electromagnetic 
modes. This and the enhancement of channeling at lower energies may help explain the 
above discrepancy. 

The possible relevance of channeling to the DAMA experiment was recently emphasized 
by Drobyshevski. ll|] He recalled the fact that ions entering single- crystals along crystal 



axes penetrate much deeper than the same ions at the same energy entering an amorphous 



form of the same material. He cited the early observations of Bredov and Okuneva[12j 
that 4 keV Cs + ions can penetrate ~ 10 3 Angstrom in Germanium crystals versus ~ 44 
Angstroms in amorphous Germanium. This "channeling" effect can significantly influence 
the ionization signals in germanium and enhance the scintillation output in Nal(Tl) single 
crystals. Specifically, channeling can dramatically increase the "quenching factor", namely, 
the ratio of the (e.m.) energy measured from a recoiling ion in the detector to that from 
an incident photon of the same energy. The quenching factor was measured in Nal(Tl) by 
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p. 0; 



a number of authors (see, eg., Ref. [13J] and most recently Ref. [141 ]) to be in the 19%-26% 
range. 

The calculation presented by Bernabei, et al.fljj] suggests that for Na or I ion recoiling 
along a channel the quenching factor may rise to be of order unity. If this is indeed correct, 
channeling which is more prevalent at lower energies generates events with higher scintil- 
lation outputs. Further, if WIMPS are light, ~ 5 Gev, then the recoil energies, ~ 1/m in 
experiments utilizing heavy nuclear targets of mass m ~Am MC i f;m are too low. However with 
the relatively light Na A ~ 26 these WIMPS can be detected in DAMA via the (channeling 
enhanced) scintillation signal, [r]] 

In this note we focus on the strong "directional" dependence of channeling. If enhanced 
channeling — accompanied by enhanced scintillation — is so important as the above might 
suggest, then the directional dependence of channeling might generate observable daily mod- 
ulations. 

More generally the prospect of a striking temporal variation of (the e.m. part of) the 
WIMP signal is interesting and may provide yet another helpful WIMP signature. 



II. The Basic Idea and Some Kinematics 

On the astronomical side our suggestion relies on the "WIMP wind" velocity W, \ W\ = 
230 Km/sec. Also the earth's rotation N-S axis makes a substantial angle of ~ 44° with W 
and the direction of W relative to the crystalline axes changes every day. 

The other physics ingredient is the direction dependence of channeling. Ions recoiling 
at angles smaller than 9(a) (or 9(p)) relative to a crystal axis or to a crystal plane are 
"Axially" or "Planarly" channeled. These "acceptance" angles decrease with increasing ion 
velocity and both planar and axial channeling probability (~ 9(p) and 9(a) 2 ) fall off as 
l/V~l/(Ef). 

For a 4 keV Iodine ion in the cubic Na-I lattice, the recipe of Appelton, et al. [17| as quoted 
in Ref. Uj yields for the main "axial" channels along the X,Y and Z directions #(a)ioo = 6.4° 
and for the other axes 9(a)uo = 4.9°; 9(a)m = 2.8°. The estimated planar channeling 
opening angles for the principle symmetry planes are: 9(p)\oo = 4.1°; 9(p)uo = 3.2° and 
9(p) in = 2.7°. 

Since the 100 planes are further apart and have a larger density of Ions the above pattern 
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is qualitatively understood. 

The above implies that an I - Ion at a random direction has a probability of: 
2tt-[2-(4.1°)x3 (for the three 100 planes) +2- (3.1°) x 6 (for the six 110 planes) +2- (2.7°) x 4 
(for the four 111 planes)] / (4tt ■ 56° ■ 4) ~ 20% 

to channel — namely, to be in the union of the cones and annular patches of all symmetry 
axes/planes. (The 56° in the denominator is a radian in degrees and the extra 1/4 is due to 
the fact that an ion incident in any direction "sees" only vr/(4 • tt) = 1/4 of the sphere.) 

We are interested not in this overall averaged channeling probability but in its directional 
dependence and the resulting temporal modulations. 

Consider Ions recoiling with a specific velocity vector V in space. As the earth rotates 
V describes a cone intersecting the unit sphere at latitude (3 = arccos[(V(||))/(|V|)] with 
V\\ the component of V parallel to the equatorial plane. As V moves in and out of any 
one of the above cones/patches, channeling along particular axes or planes turns on and off. 
Thus we expect not merely a 20% time averaged channeling probability but rather 100% 
channeling during a very specific 20% temperate subset of the day! 

To most easily visualize this, consider the artificial case (3 = 0, namely, V lies in earth's 
equatorial plane which we further take to be the crystal X-Y plane (with the 001 axis usually 
pointing along the vertical, this would happen for experiments at the north or south pole). 

Let us further assume that at 6 AM V points in the X direction. The above channeling 
acceptance angles imply that axial channeling in the 100 (+X) axis and planar channeling 
in the 010 (XZ) plane take place in the 25.6 and 16.4 minutes time intervals, respectively, 
prior to and after 6 AM. After a 45° rotation or three hours later, namely, at 9 AM +/- 
12.8 minutes, we have planar channeling in the 011 plane (at 45° with + X and +Y). At 
12 PM, the 6 AM channeling-time pattern repeats, now with the 010 Z axis and (100) YZ 
plane. At 3 PM we repeat the 9 AM pattern with the -X,Y plane, etc., so that channeling 
is modulated with a six-hour period. 

Different computable daily modulations of channeling arise for the realistic (3 ~ 45° and 
other geographical latitudes (where the artifact of having continuous channeling in the XY 
plane as in the above is avoided). We deal later with Grand Sasso's location and the actual 
(3 ~ 45° angle between W and the earth's rotation axis. 



The above is reminiscent of the Bragg enhancements at certain angles (or times) o 



makoff converted solar axions in Germanium crystals used by the Solax collaboration 18]. [19] 



Pri- 



4 



to improve the sensitivity of underground axion searches. 

Unfortunately, unlike the case where all photons point to the sun, not all ions recoil in 
the same direction! 

The velocity V of the WIMPs is the sum of wind velocity W and U, the WIMP velocity 
in the halo: 

V = W + U (1) 

For all WIMPs to be incident in the direction of W we need that (i) U be parallel to W 
and (ii) the recoil momentum q be parallel to V — the velocity of the WIMP. 

However, the halo WIMP velocity U in the halo frame is assumed to have an isotropic 
Maxwellian distribution: 

f(U)~U 2 exp-(U 2 /(U 2 )) (2) 

with Uq ~ 220 Km/sec, and also the recoil only partially tracks the direction of the WIMP. 

To find the angular distribution of the recoil ions with respect to the direction of the 
WIMP wind W we integrate over the above distribution subject to the kinematics of the 
reaction: 

WIMP (momentum p,mass M) + Ion (at rest,mass m) — > final WIMP (p') + recoil Ion(g) 
with 

p = M(W + U) (3) 

Momentum conservation then fixes p' for a given recoil q. 

The remaining relevant constraint is that of energy conservation: 

6[(p 2 -(p-qf)/(2M)-q 2 /(2m)}. (4) 

The uniform angular part of the U integration over this 5 function yields functions re- 
quiring u to be larger than the (absolute value of) 

q/n-Wcos{6) (5) 

where 9 is the angle of interest between the recoil and the WIMP wind, and fi = MmjiM + 
m) the reduced mass of the WIMP-nucleus system. Only C/'s larger than the above can 
(anti-) align with the direction of the recoil so that when added the WIND component 
W cos(#) yields the desired recoil momentum q. 
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The remaining integration over the gaussian yields the angle-momentum distribution, 

dn/(q 2 dqd cos(0)) = exp(- [<?/(2/i) - W cos(6)] 2 / (U%)) S(q 2 ) . (6) 

The factor S(q 2 ) represents nuclear form factor effects roughly parameterized by 
exp(— (q 2 ■ R 2 )) with R ~ 1.2 A 1 / 3 Fermi, the nuclear radius. For the low recoil momenta — 
dominating DAMA's events, and for which significant channeling is expected — the FF effects 
are small and can be neglected. Indeed it was claimed [3] that only relatively light ~ 5 GeV 
WIMPS with the standard halo velocity profiles are consistent with DAMA's modulation 
and all other experiments (providing that channeling is indeed important and strongly en- 
hances scintillations). In this case, \i ~ M and the average recoil momentum q = 2fiW ~ 10 
MeV with no significant FF effects, even for Iodine with R ~ 6 Fermi. The corresponding 
recoil energies are ~ 2 keV for Na and ~ 0.5 keV for Iodine. 

Using Uq ~ W since both are close to 230 km/sec, the angular distribution 

exp(-[x-cos(fl)] 2 (7) 

falls by, l/e and (1/e) 2 ~ 0.125 between 9 = and 9 = vr/2 for x = q/(2fi ■ W) = 1 and 2, 
respectively. 

For q ~ 0, W = U and there is no forward peaking of q. The overall angular distribution 
obtained by integrating Eq. ([6]) over q is reasonably peaked at 9 = 0. 

Still the width of the 9 distribution together with the high (cubic) symmetry of the Na-I 
lattice causes some planar/axial channels to be exposed to some degree to the recoiling ions 
at any time. This dilutes the dramatic modulations expected if all recoils are directed along 
the fixed W directions, but as we indicate below still leaves significant daily modulation. 

III. Possible Experimental Setups 

The latitude of Grand Sasso is ~ 43° which is close to the W direction. Hence, the WIMP 
wind hits the top 001 (XY) plane of the DAMA detectors if indeed placed vertically — almost 
head-on at the appropriate time of day on both the wide planar XZ and YZ and the X axis 
channels are then open. 

For this latitude the earth's rotation axis is close to a 101 axis, namely, the face diagonal 
of the XZ plane. Hence, unlike the simpler idealized case discussed above the lattice is not 
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FIG. 1: The angular distribution of the recoils. 

invariant under 90° rotation but rather only under a 180°. So while the system will have 
12-hour periodicity we may see now even stronger variation after six hours when we will lose 
all the above channels regaining only the 110 channel which has a smaller acceptance angle! 
The scalar product of W and W' — the original — and 90° rotated W directions is 

W-W' = 1/2 (8) 

so that the angle between them is 60°. The angular distribution of incoming WIMPs relative 
to the direction of Ws allows, even at this point in time a nonzero flux of WIMPs to fall on 
the rotated 100 plane. However the angular distribution of Fig. 1 suggests that this flux is 
reduced by a factor of 0.5. Thus 0(25%) modulations of the channeling probabilities with 
12-hour period are expected. 

Most existing or planned WIMP direct searches are at 45° ± 5° latitude so that such 
modulations (with different longitude dependent phases) are expected in all cases. 

Can modulations of channeling be translated into specific experimental signatures? 

To provide a quantitative answer we need to know how channeling of ions of given energies 
along specific axes/planes affects the intensities and/or other properties of the e.m. signal 
in any given experimental setup. 

If the channeling effects on enhanced quenching at DAMA's energies are indeed as strong 
as the analysis of Drubov and DAMA suggest — and as should be the case if channeling is 
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a key in reconciling DAMA and the other experiments — then also the direction effect and 
attendant daily variation of scintillations should be seen. 

If the individual events from the recently published DAMA/LIBRA experiment^ are 
appropriately time stamped, and the orientation of the crystals is known, the search for 
sub-daily modulation of the scintillation signal can be done with existing data! Needless 
to say, if indeed the expected daily variations with the correct phases will be seen then 
the likelihood that DAMA has indeed discovered a WIMP dark matter will be greatly 
enhanced. Using an old quotation for the DAMA data from which an annual DFS effect has 
been extracted, "There could be more gold in them mountains." 

To test the effect of different ion recoil directions the crystal can be exposed 
to a monochromatic neutron beam (with energy fixed by time of flight, or chopper 
technique). [l^], [3] From the scattered neutron the energy of the recoiling ion and its di- 
rection can be found and the procedure can be repeated for different crystal orientations. 
In thin target neutrons, like the WIMPs, collide once in the bulk. By measuring the e.m. 
response — e.g., the scintillation signal in DAMA's Nal(Tl) crystal — in each event some cali- 
bration can be directly obtained. Such calibration is indeed precise if the ratios of neutrons' 
cross sections on the two ions are the same as for the WIMPs. 

Measurements by directly bombarding with the relevant ions with the O(keV) energy are 
susceptible to strong surface effects since such ions stop in a thin < 10 3 Angstrom surface 
layer. 

We turn next to another class of detectors based on electron hole pair creation for (e.m.) 
experimental signature. In high purity germanium detectors the quenching was measured 



in the first dark matter search utilizing Ge to be ~ 25%[20|(the quenching factor here is 
defined as the ratio of charge collected in a nuclear recoil event to that induced by a photon 
of the same energy E recoi i). Usually the crystals for the high purity germanium detectors are 
pulled along the 100 axis by the Czochralski growth technique. However the remaining 110 
and 111 axes can be found by X-ray diffraction or by a technique of the XIA corporation 21] 
and the Julich group 22[ . In the experiments done at the X-ray Instrumentation Associate 
Laboratory, a weakly collimated 57 Co was used to irradiate the detector with 122 keV gamma 
rays. The source was rotated around the vertical 001 axis and the rise time of the pulses 
was meanwhile continuously measured. Maximal and minimal rise times were recorder about 
50° apart. The difference between the maximal and minimal rise times was approximately 
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20 nanoseconds, whereas the mean rise time was 190 ns. Considering the early finding of 



Bredov and Okuneva. 



121 ] channeling of the ions may strongly effect the pulse rise times and 



pulse amplitudes as well. 

The genesis of this effect is in the directional dependence of the carrier mobilities. In a 
normal p-type closed end coaxial detector, the electric fields point radially inward. Accord- 
ingly an ion liberated by a WIMP and channeled along 110 or 111 axis (i.e., not parallel to 
the cylinder axis) will travel further before stopping conceivably changing significantly the 
rise time which will thus be sensitive to the initial direction of the WIMP. 

Clearly in order to exploit this phenomenon in future Gerda and Major ana arrays should 
use one of the above X-ray or other techniques to determine the 110 and 111 direction of 
each individual crystal relative to its mounting. It would also be very valuable to perform 
the neutron scattering experiments described above and compare the pulse height and rise 
times from germanium Ions scatters in channels and the rest. References 2l|,Jul provide 
only information on the directional dependence of the carrier mobility which is related to, 
but may be quite different from, the directional effects of the WIMP scattering — though 
channeling effect for WIMPs may well be siginificant. 



IV. Summary and Conclusions 

While the present work was largely motivated by DAMA's annually varying signal and 
the potentially crucial role of Ion channeling therein, its scope is much wider. We focus on 
the strong directional dependence of channeling and its consequent variation as the crystal 
which is fixed in a terrestrial laboratory rotates relative to the fixed WIMP Wind daily. This 
in turn will modulate daily the e.m. signal since channeling appears to strongly affect the 
output e.m. signals from the detectors. It does so by modifying quenching of scintillation 
ionic crystal detectors, changing the rise time of the ionization signal in Germanium and 
potentially in many other ways in other setups. This, in turn, can lead to observed daily 
modulations in underground detectors due to the WIMP signal 

The purely kinematic DFS annual modulation is expected universally in all detectors and 
for all recoil energies. This is not the case for our suggested daily variations which require 
crystalline detectors and substantial channeling. Still our effect may complement the DFS 
modulations. 
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It is therefore highly worthwhile to look for such a signal in all crystalline detectors and, 
in particular, in DAMA where an annual signal has been seen and using the very same set 
of data for yet an independent WIMP test is naturally called for. 

We conclude by listing various, rather obvious points and finally making a few rather 
speculative additional remarks. 

(i) We focused here on low energy recoils expected for low WIMP masses. The more 
conventional LSP WIMPs have relatively high masses >~ 100 Gev, and for detector 
nuclei of similar mass we expect high, ~ 10-100 keV, generic recoil energies. 

This diminishes our effect in two ways. First, channeling is reduced at the higher 
energies. This may be mitigated to some extent by the FF effects strongly damping 
the high recoils which q « 2/j,W implies. However, in this case the x parameter in 
Eq. ([7]) above approaches zero and there is no peaking at 9 = i.e., no tendency of 
the recoils to be in the direction of wind W — further undermining our arguments. 

(ii) The high symmetry of the cubic Na-I crystal helps overall channeling by providing 
many channeling axes/planes. It does however limit the magnitude of the modulation, 
suggested as we have several open channels in the angularly broad beam of incident 
WIMPs at all times. One could conceive of using just for this purpose lower sym- 
metry lattices such as graphite with a pronounced preferred single channeling plane. 
Unfortunately most low symmetry crystals consist of "bad" detector materials and, in 
particular, are not transparent to light. 

(iii) Already some time ago, P. Sikivie suggested a complex filament-like pattern for the 
phase-space of cold dark matter. Recent large N body simulations and observations of 
Sagitarius flows suggest a complex hierarchial pattern of WIMP density and velocity 
distribution and it is possible that the local WIMP velocities differ somewhat from 
the standard scenario underlying the DFS suggestion and the putative DAMA effect. 
Still, it is most likely that the bulk of the DM is not subclustered and while some 
modification will occur the basic pattern of WIMP velocities underlying DFS and our 
suggestion is largely maintained. 

(iv) The dopant Thalium ions tend to scatter the channeling ions in the Na-I crystal and 
diminish channeling. The latter would thus be much more prevalent in a crystal of 
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high purity. The following "Wiggler" mechanism may still allow some e.m. output in 
this case. The Ion channeling with a speed of ~ 100-300 km/sec encounters crystal 
ions in the wall every ~ a ~ 2-3 Angstrom. This in turn subjects it to a time varying 
electric field of frequency / = v/a ~ 10 15 sec -1 equivalent to 4 eV. If the channeling 
path is long enough, we have some laser-like coherence and we will be able to keep 
exciting and de-exciting by spontaneous emission. (Though because of atomic energy 
scales, this holds only for the hallide ion which here is quite heavy.) 

(v) It is conceivable that channeling will reflect not only in the e.m. signal but in bolometry 
as well. The recoil speed ~ 100-300 km/sec exceeds by ~ 100 that of sound in the 
material. We will have therefore a Cherenkov-like cone of supersonic shock with a 
small opening angle, ~ c/V ~ 1/100. The sound and heat will then arrive first at one 
particular crystal face, namely, the one facing the favored channeling directions at this 
time. 
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